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CUT-OFF, CRITICAL POWER LAW, AND METASTABILITY 



DAVID A. LEVIN, MALWINA J. LUCZAK, AND YUVAL PERES 

o 

Abstract. We study the Glauber dynamics for the Ising model on the 
04 complete graph, also known as the Curie-Weiss Model. For p < 1 , we 

Q prove that the dynamics exhibits a cut-off: the distance to stationarity 

D drops from near 1 to near in a window of order n centered at [2(1 - 

P)Y^n\ogn. For = 1, we prove that the mixing time is of order n^^^. 
For yS > 1 , we study metastability. In particular, we show that the Glauber 
^-H dynamics restricted to states of non-negative magnetization has mixing 

time O(nlogn). 

d 1. Introduction 

L_j 1.1. Ising model and Glauber dynamics. Let G = (V, £) be a finite graph. 

Elements of the state space Cl := will be called configurations, and 

j>. for cr € Q, the value cr(v) will be called the spin at v. The nearest-neighbor 

O energy H(cr) of a configuration cr g {-1, 1}^ is defined by 

^ H(o-) J{v,w)o-(v)o-(w), (1.1) 



v.weV, 



f-^ where w ~ v means that {w, v] G £. The parameters J(v, w) measure the in- 

^ teraction strength between vertices; we will always take J(v, w) = J, where 

^ ^ / is a positive constant. 

^ For > 0, the Ising model on the graph G with parameter is the proba- 

^ bility measure /uonQ. given by 

g-mo-) 

where Z(fi) = J^^-^ci e'^"'-"'^ is a normalizing constant. 
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The parameter /? is interpreted physically as the inverse of temperature, 
and measures the influence of the energy function H on the probability dis- 
tribution. At infinite temperature, corresponding to yS = 0, the measure /i is 
uniform over Q. and the random variables {cr(v)}vey are independent. 

The (single- site) Glauber dynamics for p. is the Markov chain on O. with 
transitions as follows: When at cr, a vertex v is chosen uniformly at random 
from y, and a new configuration is generated from jj. conditioned on the set 

{77 6 Q : ri{w) = cr(v), w 4^ v}. 

In other words, if vertex v is selected, the new configuration will agree with 
cr everywhere except possibly at v, and at v the spin is -i-l with probability 

P^^-^^ •= (1-3) 

where S*'((t) := / Sh iw-v Evidently, the distribution of the new spin 

at V depends only on the current spins at the neighbors of v. It is easily seen 
that (Xf) is reversible with respect to the measure ju in (1.2). 

In what follows, the Glauber dynamics will be denoted by (^f)^o- 
use Po- and Eo- respectively to denote the underlying probability measure 
and associated expectation operator when Xq = cr. 

A coupling of the Glauber dynamics with starting states cr and cr is a 
process (X,, X,),>o such that (X,) is a version of the Glauber dynamics with 
starting state cr and (X,) is a version of the Glauber dynamics with starting 
state a. If a coupling (X,, is a Markov chain, we call it a Markovian 
coupling. We write Pg-o- and E^^- for the probability measure and associated 
expectation respectively corresponding to a coupling with initial states cr 
and a. 

1.2. Order n log n mixing and cut-off. Given a sequence G„ = E„) of 
graphs, we write /j.,, for the Ising measure and {X") for the Glauber dynamics 
on G„. The worst-case distance to stationarity of the Glauber dynamics 
chain after t steps is 

dn(t) := max ||P^(Xf e •) - Mnhv, (1-4) 

where \\fi - vHxv denotes the total variation distance between the probability 
measures fi and v. The mixing time ?mix(^) is defined as 

W(n) := min{f : 4(0 < 1/4}. (1.5) 

Note that ?mix(«) is finite for each fixed n since, by the convergence theorem 
for ergodic Markov chains, J„(0 ^ as ? ^ 00. Nevertheless, t^a{n) will 
in general tend to infinity with n. Our concern here is with the growth rate 
of the sequence ?mix(")- 



GLAUBER DYNAMICS IN MEAN-FIELD 



3 



The Glauber dynamics is said to exhibit a cut-off at with window {w„} 
if Wn = o(tn) and 



The first part of this paper is motivated by the following conjecture, due 
to the third author: 

Conjecture 1. Let (Gn) be a sequence of transitive graphs. If the Glauber 
dynamics on Gn has t^a{n) = 0(n log n), then there is a cut-off. 

We establish this conjecture in the special case when G„ is the complete 
graph on n vertices and /3 < 1 (the "high temperature" regime), where the 
Glauber dynamics has 0{n log n) mixing time. 

1.3. Results. Here we take G„ to be Kn, the complete graph on n vertices. 
That is, the vertex set is V,, = {1,2,..., n}, and the edge set contains all 



(^2] pairs {/, j} for 1 < / < j < n. We take the interaction parameter J to be 
1 /n; in this case, the Ising measure on {-1, 1}" is given by 



In the physics literature, this is usually referred to as the Curie-Weiss model. 
For the remainder of this paper, Ising model will always refer to the measure 
ju in (1 .6), and Glauber dynamics will always refer to the one corresponding 
to this measure. We will often omit the explicit dependence on n in our 
notation. 

It is a consequence of the Dobrushin-Shlosman uniqueness criterion that 
tmix(n) = 0(n\ogn) when yS < 1 (Aizenman and HoUey, 1987). See also 
Bubley and Dyer (1997). Our first result is that there is a cut-off phenome- 
non in this regime: 

Theorem 1. Suppose that j3 < 1. The Glauber dynamics for the Ising model 
on Kn has a cut-off at = [2(1 - fi)]~^n\ogn with window size n. 

Remark 1. Most examples of Markov chains for which the cut-off phenom- 
enon has been proved tend to have ample symmetry, for example, random 
walks on groups. Part of the interest in Theorem 1 is that the chain stud- 
ied here is not of this type, and our methods are strictly probabilistic - in 
particular, based on coupling. Recently, Diaconis and Saloff-Coste (2006) 
gave a sharp criterion for cut-off (for separation distance) for birth-and- 
death chains. 



lim lim inf dn(tn - yw„) = 1 , 



lim lim sup dn{tn + JWn) = 0. 




(1.6) 
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In the critical case yS = 1 , we prove that the mixing time of the Glauber 
dynamics is order n^^^. 

Theorem 2. If (3 = 1, then there are constants Ci, C2 > such that for the 
Glauber dynamics for the Ising model on K„, 

Cin^'^ < < C2n^'^. 

Finally, we consider the low-temperature case corresponding to /3 > 1 . 
To state our result, it is necessary to mention here the normalized magneti- 
zation, the function S defined on configurations cr by 5 (cr) := n"' ^(')- 
Also, we define the set of states with non-negative magnetization, 

D.-" ■={ojeX : S{cr)>0}. 

By using the Cheeger inequality with estimates on the stationary distribu- 
tion of the magnetization, the mixing time is seen to be at least exponential 
in n - slow mixing indeed. Arguments for exponentially slow mixing in 
the high temperature regime go back at least to Griffiths, Weng and Langer 
(1966). 

In contrast, we prove that the mixing time is of the order n log n if the 
chain is restricted to the set O.^ . To be precise, the restricted dynamics 
evolve as follows on Q^: Generate a candidate move 77 according to the 
usual Glauber dynamics. If 5 (77) > 0, accept 77 as the new state, while if 
S {rf) < 0, move instead to -77. 

Theorem 3. If [3 > 1 then there exist constants C^ijS), C^ifi) > depending 
on fi such that, for the restricted Glauber dynamics for the Ising model on 

Kn, 

C3(fi)n\ogn < t^M < C4(J3)n\ogn. 

For other work on the metastability of related models, see Bovier, Eck- 
hoff, Gayrard, and Klein (2001, 2002), and Bovier and Manzo (2002). 

The rest of the paper is organized as follows: Section 2 contains some 
preliminary lemmas required in our proofs. Theorems 1, 2 and 3 are proved 
in Sections 3, 4, and 5, respectively. Section 6 contains some conjectures 
and open problems. 

2. Preliminaries 

2.1. Glauber dynamics for Ising on K„. We introduce here some nota- 
tion specific to our setting of the Glauber dynamics for the Ising model on 
Kn. For a configuration cr, recall that the normalized magnetization S(cr) is 
defined as 

1 " 

S{cr):=-ya(j). 
n ^ 
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Given that the current state of the chain is a and a site i has been selected 
for updating, the probability pier, i) of updating to a positive spin, displayed 
in (1.3), is in this case p+(5(cr) - n~V(/)), where />+ is the function given 
by 

l+tanhpe^) 

^= ;^77^ = ^ • (2.1a) 

Similarly, the probability of updating site i to a negative spin is p-{S{cr) - 
n~^cr{i)), where 

e-^-' l-tanh065) 
^= ^^77^ = 2 • (2.1b) 

2.2. Monotone coupling. We now describe a process called the gran J cow- 
pling, a Markov chain ({X,°"}o-en)f>o such that for each cr e Q., the coordinate 
process (X'j^)t>o is a version of the Glauber dynamics started at cr. It will 
suffice to describe one step of the dynamics. Let / be drawn uniformly from 
the sites {1,2, ... ,n}, and let U he a uniform random variable on [0, 1], in- 
dependent of /. For each cr e Q, let U determine the spin S"' according 
to 

^^^1+1 0<U<p^(S((T)-n-'cTim, 
[-1 p+{S{(r) - n-^o-(I)) <U <1. 
For each cr, generate the next state X'[ according to 

zr(o = 

We write P# and for the probability measure and expectation operator 
on the measure space where the grand coupling is defined. 

For a given pair of configurations, cr and cr, the two-dimensional projec- 
tion of the grand couphng, (Xf, X'j^)t>o, will be called the monotone coupling 
with starting states cr and a. 

For two configurations cr and cr', the Hamming distance between cr and 
cr' is the number of sites where the two configurations disagree, that is 

1 " 

dist((r, (r'):=-Y, k(0 - (t\i)\. (2.2) 

^ ,-=1 

Proposition 2.1. The monotone coupling {Xt,X,) of the Glauber dynamics 
started from cr and & satisfies 

E,, [dist(X„ 1,)] < p'dist(cr, &), (2.3) 

where 

p:=l-n-\\-n tanh(;0/n)) . (2.4) 
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Proof. We first show that (2.3) holds with t = \ provided dist(cr, cr) = 1. 
Indeed, suppose that cr and & agree everywhere except at /, where a{i) = - 1 
and (t(/) = + 1 . 

Recall that the vertex which is updated in all configurations in the grand 
coupling is denoted by /. If / = /, then the distance decreases hy \ \ ii I i 
and the event 5(7) occurs, where 

BU) := {pAS{a) - o-(j)/n) <U< p^(S{d-) - a(j)/n)} , 

then the distance increases by 1 . In all other cases, the distance remains the 
same. Consequently, 

dist(Xi,li) = 1 - 1{/ = /} + 1{/ = (2.5) 

j*' 

Note that 5((t) - d-(j)/n = 5(cr) - o-(j)/n + 2/n for j 4^ i. Thus, letting 
Sj/n = 5(cr) - o-U)/n, for j 

P^(5(j)) = ^ [tanh(J3(Sj + 2)/n)) - tanh(J3sj/n)] < tanh(J3/n). (2.6) 

Taking expectation in (2.5), by the independence of U and / together with 
(2.6), 

Eo.[dist(Zi,Xi)] < 1 - - + tanh(/3/n) = p (2.7) 

n 

This establishes (2.3) for the case where cr and cr' are at unit distance. 

Now take any two configurations cr, (t with dist(cr, cr) = k. There is a se- 
quence of states cro, . . . ,crk such that ctq = cr,cri, = a, and each neighboring 
pair cr,, cr,_i are at unit distance. Since we proved the contraction holds for 
configurations at unit distance, 

k 

[dist(Zf ,Xf)] < [dist(Xf',X;^'-')] <pk = pdist((r, a). 

i=\ 

This establishes (2.3) for ? = 1; iterating completes the proof. ■ 

We mention another property of the monotone coupling, from which it 
receives its name. We write a < a' io mean that cr{i) < cr'{i) for all i. 
Given the monotone coupling (Xt,Xt), if Xt < X,, then X^ < Xs for all s > t. 
This is obvious from the definition of the grand coupling, since the function 
p+ is non-decreasing. 

2.3. Magnetization chain. Let 5, := S(Xt), and note that (5,) is itself a 
Markov chain on Q5 := {-1, -1 -I- 2/n, . . . , 1 - 2/n, 1}. The increments 
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St+i -St take values in {-2/n, 0, 2/n}, and the transition probabilities are 



Puis, s') 



^^p-(s-n-^) s' = s-2/n, 

+ s' = s + 2/n, (2.8) 

I- ^^p^is-n-^)- ^p^is + n'^) s' = s, 

for 5 6 Q5, where p+{s) and p^{s) are as in (2.1). 

Remark 2. It is easily verified that Pm{-s, -s') = Pm(s, s'), so the distribu- 
tion of the chain (5 ,) started from s is the same as the distribution of (-5 ,) 
started from -s. 

Remark 3. Let (X,^) be the Glauber dynamics restricted to Q^, and define 
S'l := 5(X/'). The chain (5/^) has the same transition probabilities as the 
chain 

In the remainder of this subsection, we collect some facts about the 
Markov chain (St) which will be needed in our proofs. 

If (Xt, Xt) is a coupling of the Glauber dynamics, we will always write St 
and St for S{Xt) and 5(X,), respectively. 

Lemma 2.2. Let p be as defined in (2.4). If{Xt,Xt) is the monotone cou- 
pling, started from states a and cr, then 

E^,^ [|5, - 5,1] < |^jp'dist((r,a-) < 2p\ (2.9) 

Proof. Using the triangle inequality, we see that jS'f-S,! < (2/n)dist(Xf, X^. 
An application of Proposition 2.1 completes the proof. ■ 

Lemma 2.3. For the magnetization chain (S ,), for any two states s and s in 
Qs with s > s, 

0<E,[5i]-Ej[5i] <p(s-s). (2.10) 

Also, for any two states s and s, 

|E,[5i]-E,,[5i]|<p|5-S|. (2.11) 

Proof Let (X,,X,) be the monotone coupling, started from (cr, (t), where 
cr > a- and 5(cr) = s,S(d-) = s. In this case, s - s = (2/n)dist(cr, (t), and 

^crAlSi -Si\] = E,,^[(2/«)dist(Xi,li)] < -pdist((r,a-) = p(s - s). 

n 

By monotonicity, Xi > Xi and so 5i > 5i. Thus, Eo-Li*!] - Ea-[5i] = 
Eo-,a-[|5' 1 - i|], which, together with the preceding inequality, proves that 

EJSi]-E^[Si]<p(5-5). (2.12) 
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The left-hand side of (2. 12) equals £^[5 i]-Ej[5 1], because (5 ,) is a Markov 
chain. Moreover, the left-hand side does not depends at all on the coupling. 
This proves (2.10). An analogous bound in the case S((t) > 5(cr) estab- 
lishes (2.1 1). ■ 

We now study the drift of (St) in some detail. From (2.8), 

2(l-s\ , 2/l + 5\ , 

E[Sr^,-Sr\Sr = s] = - p^{s + n') - - \—-\ p^{s - n'), 

and hence 

E[5,^i - 5, I 5, = ^] = - [Us) - s + Onis)] , (2.13) 

n 



where 



Us) := ^ {tanh[^(5 + n')] + ismfi{s - n')]] 

0^(s) := y {tanh[y3(5 + n"^)] - tanh[J3(s - n"')]) . 
The approximation 

E[5,+i - 5, I 5, = 5] « - [tanhOe*) - s] (2.14) 

n 

will play an important role in our proofs, and we will need to control the 
error fairly precisely. For the moment, let us observe that (2.14) is valid 
exactly as an inequality for s > 0: 

E[St+i-St\St = s]<- [tanhCe*) - s] . (2.15) 

n 

This follows from the concavity of the hyperbolic tangent, together with the 
fact that the term Onis) in (2.13) is negative. By Remark 2, for < 0, 

E[5^+i - 5^ I = 5] > - [tanhCe^) - s] . (2.16) 

n 

Since (St) does not change sign when > n~\ and because tanh is an odd 
function, putting together (2.15) and (2.16) shows that, for > n~^, 

e[|5,+i| I 5,1 < 15,1 + - [tanh(^|5,|) - |5,|] . (2.17) 

Since tanh(.x;) < x for x > 0, when /3 < I equation (2.15) implies that, for 
s>0, 

E[5,^i-5,|5, = ^] < (2.18) 

n 

Define 

To := mf{t>0 : |5J < l/n}. (2.19) 
Note that, for n even, |5rJ = 0, while for n odd, |5rol = l/«- The notation 
To will be used with the same meaning throughout the paper. 
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On several occasions, we will need an upper bound on the probability 
that an unbiased random walk remains positive for at least u steps. The 
following lemma gives a classical estimate. 

Lemma 2.4. Let {Wt)t>o be a random walk with E[V7,+i - Wt \ V7,] = and 
\Wt+\ — W,\ < B for some constant B. Then there is a constant c > such 
that, for all u, 

P.dW^il >0,...,|WJ >0)< ^. (2.20) 

(Here indicates probabilities for the random walk started with Wq = k.) 

Lemma 2.4 can be proved using hitting estimates in Feller (1971); al- 
ternatively, it can be seen to be a special case of equation (3.9) in Bender, 
Lawler, Pemantle, and Wilf (2003/04). 

The following lemma is proved for n even. The proof can be modified to 
deal with the case of n odd by replacing with \/n;we omit the details. 

Lemma 2.5. Let fi < I, and suppose that n is even. There exists a constant 
c such that, for all s and for all u,t >0, 

cn\s\ 

P( 15 J > 0, . . . , > I 5„ = ^) < (2.21) 

Proof. It will suffice to prove (2.21) for 5 > 0, in which case the absolute 
values may be removed. 

By (2.18), E[5,+i - 5, | 5,] < for 5, > 0. Also, there exists a constant 
b > such that P(5;+i - Sti'0\St)>b for all times t, uniformly in n. It 
follows that (St) can be coupled with an unbiased nearest-neighbor random 
walk (Wt) on Z satisfying 

• F(Wi - ^ I W^o = w) = Z7 for all w, 

• Wo = ns/2, 

• nSt/2 < Wt for t less than the first time u when 5„ < n~^. 
From Lemma 2.4, there exists a constant c > such that 

P(5„+i > 0, . . . >0\Su = s)< P(Wi > 0, . . . , W > I Wo = ns/2) 

ens 
< . 

■ 

2.4. Variance bound. 

Lemma 2.6. Let (Z;) be a Markov chain taking values in R and with tran- 
sition matrix P. We will write P^ and for its probability measure and ex- 
pectation, respectively, when Zq = z. Suppose that there is some < p < 1 
such that for all pairs of starting states (z, z), 

|E,[Z,]-E,[Z,]|<pqz-2|. (2.22) 
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Then Vt := sup^^ Var^o(Zf) satisfies 

Vt < vimin{f,(l - p^Y^}. 

Remark 4. Suppose that, for every pair (z, 2), there is a coupling (Zi, Zi) of 
P{z, •) and ^(2, •) such that 

E,4|Zi-Zi|]<p|z-2|. (2.23) 

By iterating (2.23), 

|E,[Z,] - Ej[Z,]| < E,,,[|Z, - Z,|] < p'\z - z\. 

The left-hand side does not depend at all on the coupling, and in particular, 
(2.22) holds. Moreover, if the state-space of (Z,) is discrete with a path 
metric and (2.23) holds for all neighboring pairs z,z, then it holds for all 
pairs of states; see Bubley and Dyer (1997). 

Proof. Let (Z,) and (Z*) be independent copies of the chain, both started 
from Zo- By the Markov property and (2.22), 

I E,„[Z, I Zi = zi] - E,„[Z; I Z* =z\]\ = \ E„ [Z,_i] - E,* [Z* J | 

<p'-\zi-z\\. 

Hence, letting ip{z) = E,[Z,_i], we see that 

Var,„ (E,„[Z, | Z,]) = ^E,„ [[^(ZO - ip{Zt)f] 

< vyP^^'-^\ (2.24) 
By the "total variance" formula, for every zo, 

Var,„(Z,) = E,„ [Var,„(Z, | Z,)] + Var,„ (E,„[Z, | Z,]) , 

so that 

V, < sup {E,„[Var,„(Z, | Z,)] + Var,„ (E,„[Z, | ZJ)} . (2.25) 

zo 

Now, Varj.Q(Z, | Zi = zi) < v,-i for every zi, and so 

E,jVar,„(Z,|Zi)] <v,_i. (2.26) 
Thus we have shown that v, < Vt-\ + vip^'''~^\ whence 



Vt < vi Y.P'^'''^ ^ ^1 min I (1 - p^rK t] . 



i=0 



Proposition 2.7. Iffi < I, then Var(5f) = 0(n ^) as n ^ oo. Iffi = 1, then 
Var(5,) = 0(t/n'^) as n ^ ^. 
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Proof. The conclusion follows from combining Lemma 2.3 with Lemma 
2.6, and observing that vi is bounded by {A/nf since the incrememnts of 
{St) are at most 2/n in absolute value. ■ 

2.5. Expected spin value. In order to establish the cutoff at high tempera- 
ture, not only do we need to consider the magnetization chain, but also the 
number of positive and negative spins among subsets of the vertices. 

Lemma 2.8. Let [5 < 1. 

(i) For all cr e Q and every i = 1,2, ... ,n, 

\EASt]\ < le-'^-f^^"", and |E^[X,0-)]| < le-'^-f^^'^" . 

(ii) For any subset A of vertices, if 

M,(A):=^-yX,(i), (2.27) 

then \Eo-[M,(A)]\ < lAle'*-^'^^*'" andVaiiMfiA)) < cn for some constant 
c> 0. 

(iii) For any subset A of vertices and all cr e CI, 

[\MriA)\] < ne'^^-f^^"" + 0{ V")- (2.28) 

Proof. Let 1 denote the configuration of all plus spins, and let {Xj ,X,) be 
the monotone coupling with = 1 and such that Xq has distribution /i. 
(Note that then X, has distribution /i for all ? > 0, by stationarity.) From 
Lemma 2.2, because E^[5f] = 0, we have 

El < El,, [\S] - 5,1] + E, [S] < 2e-^^'-P^I\ 
By symmetry, Ei (O] ^ 2e~^^~^^'^" for all i. By monotonicity, for any cr. 

Because the chain {-S ,) started from -cr has the same distribution as the 
chain (5 ,) started from cr, 

_2e(i-W« <E^[X,(0]- 

For part (ii), the bound on the expectation follows from (i). As for the 
variance, since the spins are positively correlated, 

/ n 



Var 



V ieA 



y.xm <var Yxm 



V ;=1 



<n2Var(5,)<cn, (2.29) 



where the last inequality follows from Proposition 2.7. 

For part (iii), let {Xt,Xt) be the monotone couphng with Zq = cr and the 
distribution of Xq equal to /i. From the triangle inequality, 

E.[|M,(A)|] < E^,, [\M,{A) - M,{A)\] + E, [|M,(A)|] . 
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By the Cauchy-Schwartz inequality and since |M((A) - A/,(A)| < dist(X,, X,), 

E..[|M,(A)|] < E^,^ [dist(X„l,)] + ^E^ [MtiAf]. 
Applying Proposition 2.1 shows that 

E.[|M,(A)|] < np' + ^E,[Mr(Ay]. (2.30) 
Since the variables {X,(z)}"^j are positively correlated under fi, 

2 2 

[m,(A)2] < ^E^ [St] = j Var^(S,) = 0{n), (2.31) 

where the last equality follows from Proposition 2.7. Using (2.31) in (2.30) 
shows that 

E^ [\Mr(A)\] < ne-^^-^^"" + 0( yfn). (2.32) 

■ 

2.6. Coupling of chains with the same magnetization. The following 
lemma holds at all temperatures, though we will only be using it for /3 > 1. 
It shows that once the magnetizations of two copies of the Glauber dy- 
namics agree, the two copies can be coupled in such a way that the entire 
configurations agree after at most another 0{n\ogn) steps. Note that this 
simple coupling is not fast enough to show cutoff (where we need that once 
the magnetizations agree, only order n steps are required to fully couple). 
A more sophisticated coupling for this purpose is given in Section 3. 

For any coupling {X,,X,), we will let r denote the coupling time: 

T := mini? > : X, = 1^}. 

Lemma 2.9. Let (T,d- e Q be such that S (cr) = S (cr). There exists a cou- 
pling (Xf, Xf) of the Glauber dynamics with initial states Xq = cr andX^ = & 
such that 

lim sup Po-,a-(T > co(J3)n log n) = 0, 
for some constant co(fi) large enough. 

Proof. To update the configuration X, at time t, proceed as follows: Pick a 
site / 6 {1, 2, . . . ,n} uniformly at random, and generate a random spin S 
according to 

^ _ 1 + 1 with probability p+{S t - Xf{I)/n), 
1-1 with probability p-{S t - Xt{I)/n). 

Set 

16 1 = 1. 
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As for updating Xt, if Xt{I) = Xtil), then let 



\Xr(i) i + I, 
\S i = I. 



If Xt{I) 4^ X,{I), then we pick a vertex / uniformly at random from the set 
{/ : 1,(0 ^ Ui), and 1,(0 = X,(I)}, 

and set 



\S i = I. 



Let Dt = 2iLi \Xtii) - Xt{i)\/2 be the number of differing coordinates be- 
tween Xt and X,. 

There exists a constant ci = ci(J3) > such that p+(s) A p-(s) > ci 
uniformly over all 5 6 { - 1 , . . . , 1 } and all n. If X,(/) = 1,(7), then -D, = 
while if X,(/) ^ 1,(7), then A+i - 7), = -2. It follows that 

2ciA 

E[A+i-A|X„X,]< 

n 

so = Dt{l - Icjny' is a non-negative supermartingale, whence 



E[A] < E[7)o] 1 1 - ^ I < ne 



Taking f = conlogn for a sufficiently large constant cq = c^ifi), we can 
make the right hand side less than 1 /n, say. Markov's inequality yields 

Po-(t > Co/1 log n) < Po-(Dconlogn ^ l) ^ Eo-[7),.o„iog„] < -. 



3. Cutoff for the Glauber dynamics at high temperature 

In this section we prove Theorem 1. As always, (Xt) will denote the 
Glauber dynamics, and St = S(X,) = Yli=\ Xt{i) is the normalized mag- 
netization chain. Recall the definitions 

tn = [2(1 -l5)r'n\ogn, 
p=\-{l-l5)ln. 
To = mm{t > : |S,I < 1/n}. 
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3.1. Upper bound. For convenience, we restate the upper bound part of 
Theorem 1: 

Theorem 3.1. IffS < 1, then 

lim limsup <i„ ([2(1 -/3)]"^ nlogn + yn) = 0. (3.1) 

Our strategy is to first construct a coupling of the dynamics so that the 
magnetizations agree with high probability after + 0{n) steps. 

Lemma 3.2. Let cr and & be any two configurations. There is a coupling 
{X,, Xf) of the Glauber dynamics with Xq = cr and Xq = a such that, if 

Tmag := min{? > : 5, = 5,}, (3.2) 

then for some constant c > not depending on cr, & or n, 

c 

Pcr,(T(Tmag > t,, + jn) < —. (3.3) 

Vr 

Proof. Assume without loss of generality that S{cr) > 5(cr). Let (Xt,X,) 
be the monotone coupling of Section 2.2. Define := (n/2)|5, - By 
Lemma 2.2, for some ci > 0, 

E^,* [A J < ci V^. (3.4) 

Define ti := min{t > tn : |A,| < 1}. For tn < t < ti, allow {Xt) and {X,) to 
run independently. 

Since St>S,fovt<Ti, from Lemma 2.3, the process (S, - 5f),„<r<r, has 
non-positive drift. Moreover, since (Xr)f„<f<T, and {X,)t„<t<Ti are independent 
given for t > t„ the conditional probability that 5, - 5, is non- 

zero is bounded away from zero uniformly. Thus there is a random walk 
{Wt)t>t„ defined on the same probability space as {Xt,Xt) and satisfying: the 
increments 1/7;+ 1-V7, are mean-zero and bounded, n(5f- 5,) < W,on\tn,Ti), 
and n{S , - 5 , ) = W^, . 

By Lemma 2.4, 

P^,^(ri >t„ + rn\ X,„lj < P^AW,M > 0, . . . , Wr„^yn > | X,„,lj 

^ n\St,-StJ 

Taking expectation above, (3.4) shows that 

PcrAri >tn + yn)<0(r~"^). 

The number of plus spins in Xr^ is either one more than, or the same as, the 
number of plus spins in X^.^ . Match each plus spin in Xr, with a plus spin in 
Xt^ , and match the remaining spins arbitrarily. From time ti onwards, run 
a modified version of the monotone coupling, where matched vertices are 
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updated together in the two chains. Define dist' as the number of disagree- 
ments between matched vertices. The conclusion of Lemma 2.2 now holds 
for this modified monotone coupling, with the distance dist' replacing dist 
in (2.9). Thus, 

— Eo-,o-[A^|+y'„ I XtijX-Tj] 
n 

(l-y8)r' 



< e 



We conclude that 

P^.<r(rmag < + 7" + /«) > 1 - O (y"^^^) . 

■ 

3.2. Good starting states. To show the cut-off upper bound, we will start 
by running the Glauber dynamics for an initial burn-in period. This will en- 
sure that the chain is with high probability in a 'nice' configuration required 
for the coupling argument in Section 3.3. The following lemma is required: 

Lemma 3.3. For any a subset Qq c Q, 

d{tQ + = max ||P^(X,„+, e •) - ttHtv 

< max ||P,„(X, 6 •) - ttIItv + max P,(X,„ i Dq)- (3.5) 

Proof. For A c Q., we can bound \Va-{X,^+, 6 A) - n{A)\ above by 
[Pcr(X,„+, 6 A I = cro) - n{A)\ P^(X,„ = a^) 

+ [P.(X,„« 6 A I ^ Qo) - ^(A)] P.(X,„ i rio) 

Using the triangle inequality, the preceding displayed quantity is bounded 
above by 

^ |P^(X,„+, e A I = cro) - ;r(A)|P^(X,„ = (Tq) + P^(X,„ ^ Qq). 

o-oeno 

Taking a maximum over subsets A shows that 
||P^(Z,,+, 6 •) - ttIItv 

< 2 \\PAXt,^t 6 • I = (^o) - ;r|lTvP^(Xro = o-o) + P<.(X,„ ^ Qo). 
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By the Markov property, Pa-iXt^+r e ■ \ Xt^ = ctq) = Po-o(Xt e •), and 
bounding the average above by the maximum term yields 

||P^(Z,„+, 6 ■) - ^IItv < max \\P^,(Xr e ■) - n\hy + P^(X,, i Qq). 

Taking a maximum over cr e H establishes (3.5). ■ 
In the proof of Theorem 3.1, we apply Lemma 3.3 with 

Qo = {(^6 : \S{cr)\ < 1/2}. 
For a configuration ctq e Q. define 

uo := \{i : cro(z) = 1}|, vo := \{i ■ o-o(i) = 

the number of positive and negative spins, respectively, in ctq. Also, define 
Ao := {(m, v) : n/4<u,v < 3n/4}. Note that 

ctq e Qo if and only if (uq, Vq) £ Aq. (3.6) 

By Lemma 2.8, there is a constant > such that |Eo-[56Iq„]| < 1/4, 
whence, for n large enough, 

PAXeoni^o) = PAIS eon\> 1/2) 

<P^(\Se,n-EASeon]\ > 1/4) 
< \6Yar AS eon) = Oin-'). (3.7) 
The last equality follows from Proposition 2.7. 

3.3. Two-coordinate chain. Fix a configuration cfq e Q.q. For cr e Q., 
define 

U^,(cr) := |{/ e {0, 1, . . . ,n} : crii) = crod) = 1}| 
V^,{(T) := |{/ e {0, 1, . . . ,n} : cr(i) = crod) = -1}|. 

In what follows, we shall usually omit the subscript, writing simply U((t) 
for U^,((t) and V{cr) for y^„(cr). 

For a copy of the Glauber dynamics (X;), the process ([/,, Vt)t>Q defined 

by 

Ut = UiXt), and = V(Xt) (3.8) 

is a Markov chain on {0, 1, . . . , mq} x {0, 1, . . . , vq} (with transition proba- 
bilities depending on the designated configuration ctq). We will refer to 
the chain (Ut, Vt) as the two-coordinate chain, and its stationary measure 
will be denoted by n2. Note also that (Ut, Vt) determines the magnetization 
chain, as we can write 

^^^2(Ut-Vt)_u^_ (3.9) 
n n 
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It turns out that, by symmetry, the distance of the law of Xt to ju equals the 
distance of the law of (Ut, Vt) to n2, as established in the following lemma: 

Lemma 3.4. If(Xt) is the Glauber dynamics started from ctq and (U,, Vt) is 
the chain defined by (3.8) started from (Uq, Vq), then 

\\F^,(Xr 6 •) -juIItv = \\Fiu,,vo)i(Ut, Vt) 6 •) - ^2||tv. (3.10) 

Proof. Let 

Q(m,v) := {cr 6 Q : (U(o-),V(o-)) = (u,v)}. 
Since both //(■ | Q.(u, v)) and 

P^,(Xte-\(Ut,Vt) = (u,v)) 
are uniform over Q.{u, v), it follows that 

P,„(X, = 77)-M?7) 

Z1{t] e Q.(u, v)} r„ 
' [Vcrom, V,) = (U, V)) - fi{n(u, V))] . 

£2(m, v) 

u,v 

Applying the triangle inequality, summing over rj, and changing the order 
of summations shows that 

||P^„(X, 6 •) - juIItv < l|P(r,o,vo)((C/., V,) e •) - ^2IItv. 
The reverse inequality holds since (Ut, Vt) is a function of (Xt). ■ 

Identity (3.10) implies that it suffices to bound from above the distance 
to stationarity of the two-coordinate chain. 

Lemma 3.5. Suppose two configuration cr and & satisfy S (cr) = S (a) and 
Rq = U(d-) - U((t) > 0. Define 

Si := [o- : min{[/(cr), uo - U(cr), V(cr), v,, - V(cr)} > n/16}. (3.11) 

There exists a Markovian coupling (Xt, Xt) of the Glauber dynamics with 
starting states Xq = cr and Xq = a such that the following hold: 

(i) S(Xt) = S(Xt)forallt>0. 

(ii) IfRt := U(Xt) - U(Xt) and Ro > 0, then Rt>0 and for all t and 

E^,^[Rt^i-Rt\Xt,Xt]<0. (3.12) 

(iii) There exists a constant c not depending on n so that on the event {Xt e 

ill, Xt e 

F^^dRt^i-Rti^O\Xt,Xt)>c. (3.13) 
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Figure 1. The vertices in Xt and X, are partitioned into four categories. 



Proof. Given the coupling (X,,Xt), we define Ut := U(Xt) and V, := V(Xt), 
and note that Ut = Ut + Rt and Vt = Vt + Rf 

For any configuration cr, we divide the vertices into four sets: 



A(o-) = {i 6 


{1,2,. 


.n} 


■ o-o(i) = 


+ 1, cr(/) = 


+1}, 


B{(r) = {i 6 


{1,2,. 


.n} 


■ cro(z) = 


+ 1, cr(0 = 


-1}, 


C((r) = {/ 6 


{1,2,. 


.n} 


: CToii) = 


-1, (T(i) = 


+1}, 


D((r) = {i 6 


{1,2,. 


.n} 


■ o-oii) = 


-1, (T(i) = 


-1}, 



(3.14) 



and so 

\A((t)\ = U{(T\ \B{(T)\ = uo- U{cr), \C{cr)\ = vo - V{(t), \D{cr)\ = V{ct). 

See Figure 1 for a schematic representation of this partition for Xt and Xt. 

Our coupling is as follows: To update X,, select a uniformly random / e 
{ 1 , 2 . . . , r}, and generate a random spin S for / according to the distribution 

^ 1 + 1 with probability p+(Sf - 

I - 1 with probability p-{S t - Xt{I)/n). 

Set 

^-^'^ = 1^ i = I. 
For Xf, we select / uniformly at random from {i : Xt(i) = Xt(I)}, and let 



16 1 = 1. 
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The difference Rt+i - Rt is determined by the values of /, / and S according 
to the following table: 



/ / 


S 


Rnl - Rt 


/ 6 B(X,) I 6 D(Xt) 


+1 


-1 


IeC{Xt) leAiXt) 


-1 


-1 


l€A(X,) l€C(Xr) 


-1 


+1 


leDiXt) 


+1 


+1 


all other combinations 








It follows that 

P^,^(i?,^i -Rt = -1 I Xt,X,) = a{Ut,Vr,Rt), 
F^ARt^, -Rt = +1 I Xt,X,) = b(Ut,Vr,Rt), 
where (using the identities Ut = Ut + Rt and Vt = Vt + Rt) 

a(Ut,Vt,Rt) = (hzli)i _ - 1/n) 

b(Ut,Vt,Rt) = {-)i'^^^^f^]p-(St-Un) 
\n / \vo + U,- Vt) 

\nl\uo-Ut + Vtj 

We obtain 

E,,,^ [r,^, -Rt I Xt,x] = b(Ut, Vt,Rt) - a(Ut, Vt,Rt) 

= —[p.{St-l/n)+p^{St+l/n)], 
n 

so, in particular, 

¥.^M^,-R,\X,X]<Q. (3.15) 
Furthermore, on the event [X, e Si,Xf e Si}, 

P^,^(i?,+i -Rt^O\ XtX) > b(Ut, Vt,Rt) > c 

for some constant c > 0, uniformly in n, since the functions p+ and p are 
uniformly bounded away from and 1 . ■ 

Proof of Theorem 3.1. Applying Lemma 3.3 with = Bon, together with 
the bound (3.7), shows that 

dniO^n + t)< max ||P^„(X, 6 •) - A/|Itv + Oin'). (3.16) 
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Hence, using Lemma 3.4 and (3.6), 

dniOon + t)< max \\F^a^j,J(U„ Vr) e ■) - n2\\jw + (3.17) 

(uo,vo)eAo 

recalling that Aq = {(w, v) : n/4 < u,v < 3n/A}. 

We will call a pair of chains (Ut, Vt)t>o and (f/,, Vt)t>o a coupling of the 
two-coordinate chain with initial values (uq, vq) and (u, v) if 

• The two chains are defined on a common probability space, 

• Each of (Ut, Vt) and (Ut, Vt) has the same transition probabilities as 
(U(Xt), V(Xt)), where (Xt) is the Glauber dynamics, 

• (Uo, Vo) = (Mo, vo) and (U, V) = (u,v). 

We will always consider couplings which have (uq, vq) 6 Aq, but (u, v) will 
not be so constrained. 

For a given coupling of the two-coordinate chain as above, we let 

T, :=minU>0 : (Ut,Vt) = (Ut,Vt)}. 

For a coupling with initial states (Uq, vq) and (u, v), 

l|Pr,o,vo ((U„ V,) 6 ■) - Kv ((U, Vt) e •) IItv < P(.o,vo),(«,v)(^. > 0. (3.18) 

(See, for example, Lindvall (2002, Equation 2.8).) A simple calculation 
shows that 

max \\Puo,vo((Ut, Vt) 6 •) - ttiIItv 

(i(o,vo)eAo 

< max |iP,„,-,((t/„y,)6-)-Ps.v((f/r,Vr)6-)llTv. (3.19) 

(Mo,vo)eAo, 
(u,v) 

We say that f(n, t) is a uniform coupling bound if for any initial states 
(mo,vo) 6 Ao and (u,v), there is a coupling of the two-coordinate chain 
with 

P(r,o,vo),(«,v)(Tc >t)< f(n,t). 
If f(n,t) is a uniform coupling bound, then combining (3.18) with (3.19) 
shows that 

_rnax ||"«o,vo 

((Ut,Vt)e-)-n2\\Tw<f(n,t), 

(no,vo)eAo 

and by (3.17), 

d„(Oon + t) < f(n, t) + 0(n~^). 

Recall that t,, = [2(1 - /3)]"'(nlogn). For any y > 0, let tn(r) := t„ + yn. 
The theorem will be proved if we can establish a uniform coupling bound 
f(n, t) such that 

lim lim sup f(n, tn(y)) = 0. 

Fix (uo,Vo) 6 Aq and arbitrary (u,v). Let ctq be any configuration with 
(U(o-o), V((To)) = (uq, Vq), and let cr be any configuration with (U (cr), V(d-)) = 



GLAUBER DYNAMICS IN MEAN-FIELD 



21 



(m, v). We will construct, in two phases, a coupling {Xt,Xt) of the full 
Glauber dynamics with initial states Xq = ctq and Xq = a. Given such a 
coupling, the projections 

([/„ Vd := (UiX,), V(X,)), and V,) := (C/(X,), V(X,)) 

are a coupling of the two-coordinate chains, started from (uq, Vq) and (u, v). 

The magnetization coupling phase, lasting from time to time ?„(y) will 
ensure that St„(j) = St„(y) with high probability, and that 



During the two-coordinate coupling phase, from time t„{y) to time ?„(2y), 
with high probability the chains (Ut) and (U,) coalesce. To facilitate coales- 
cence, we must ensure that throughout the second phase with high probabil- 
ity X, 6 S] and 6 Si, where Si is as defined in (3.1 1). Also, the coupling 
will ensure 5, = 5, for all ? 6 [f„(y), t„{2y)]. 

(i) Magnetization coupling. Recall that r^ag, defined in (3.2), is the first 
time the normalized magnetizations agree. Let Hi := {Tmag ^ tniy)} be the 
event that the magnetizations couple by time t„{y). By Lemma 3.2, there 
exists a constant c not depending on ctq or cr such that 

(ii) Two-coordinate chain coupling phase. Assume that Ut^ > Ut„', if this 
is not the case, just reverse the roles of X, and X, in what follows. On the 
event H\, for t > t„{y), use the coupling constructed in Lemma 3.5. On the 
event //j, we let the two chains run independently for t > tniy). 

The outline of the remainder of the proof is as follows: By (3.12), the drift 
of the difi'erence Ut - Ut is non-positive, so it can be dominated by a process 
with independent and unbiased increments with values in {-1,0, 1}, until 
Ut - Ut hits zero. Provided that the increments of f7, - Ut are non-zero with 
probability bounded away from uniformly in n, the dominated process 
can be taken to be an unbiased random walk. We will establish that at time 
t„{y), the beginning of the second coupling phase, the expected difi'erence 
'E'(TQ,AUt„(y) - Ut„iy)] is order y/n. Thus by comparison with random walk, 
the two-coordinate process will couple in 0{n) more steps. 

We begin by showing that, if H2{t) := {Xt e Si, e Si}, then 



p 

-T 0-0,0- 



/ 

1 



U H2{ty 

\t„(y)<t<t„(2y) 



0{n'). (3.20) 



(Note that the bound above depends on y. This does not pose a problem, 
because the limit in n is taken before the limit in y in (3.1).) 
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Recall the definition of M,(A) in (2.27). We introduce the following def- 
initions: 

Ao := {/ : cToii) = 1}, 

B*:= [j {\MM\>n/32}, 

te[t„+yn,t„+2rn] 

Y := Mm(Ao)\ > n/64}. 

te[t„+yn,t„+2yn] 

(Note that |Ao| = uq.) Since M,(Ao) has increments in {-1, 0, 1}, if |M,„(Ao)| > 
n/32, then \M,{Aq)\ > n/64 for all t in any interval of length n/64 containing 
to. Consequently, B* c {¥ > n/64-} and 

P.o,^(5*) < P^,AY > n/64) < ^2^^2:£Z1. 

n 

By Lemma 2.8(ii), P^o,^(|M,(Ao)| > n/64) = 0(n-^) for t > t,,, so E^^Y] = 
0(1) and 

P^.AB*) = 0{n-'). 
Making analogous definitions and deductions for the chain {%) shows that 

V^.AB*) = 0{n-'). 

If Ut < n/ 16, then uq - Ut > 3n/l6, since we are assuming that uq > n/4. 
Consequently, if i7, < n/16, then 

|M,(Ao)| = \Ut - iuo - U,)\ > iuo - Ut) -Ut>'^. 

Similarly, uq - Ut > n/\6 implies that |M,(A())| > 1/8. An analogous ar- 
gument applied to V, and vq - Vt shows that if either Vt or vq - Vt does not 
exceed n/16, then |M,(Ao)| > n/8, since \Vt - (vq - Vt)\ = |vol > n/4. Fi- 
nally, the same implications are obtained for the chains (Xt), (Ut) and (Vt). 
To summarize, 

H2(tY c {|M,(Ao)l > n/16} U {\MtiAo)\ > n/16}. 

Thus, 



p 

0-0,0- 



U H^it)' 



_<P^„,^(5*) + P^„,^(5*) = 0(n-i). 

\t„{y)<t<t„{2y) 



Recall that Rt = \Ut - and let H2 := r\t„{y)<t<t„{2y) ^2(0- On the event 
H2, the process Rt can be dominated by a nearest-neighbor random walk, 
with delay, until the first time when (R,) visits 0. Then by Lemma 2.4, on 



ny 
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Taking expectation gives 



P.o,^ (.{rc > tn(2r)] n //2 n H^) < (3.21) 



Cl'Eo-Q,o-[\Rt„(y)\] 

Observe that 

Ut = M,{Aq) + uq/2, and Or = M,(Ao) + uo/2, 

whence 

\Ur - Ur\ = \M,(Ao) - MtiAo)\ < \Mr(Ao)\ + \Mt(Ao)\. 

Taking expectation shows that 

E.o,*[l^.l] < E,„[|M,(Ao)|] +E^[|M,(Ao)|]. 

Applying Lemma 2.8(iii) shows that Eo-(,,o-[|i?,„(y)|] = 0( V") • 
Using this estimate in (3.21), we conclude that 

r 0-0,0- 

(r, > tnily)) < ({r, > t„(2y)} n //2 n //i) 
< ^ + 0{n-'). 

This gives the uniform coupling bound required. ■ 

3.4. Lower bound. Recall = [2(1 - /3y^]n log n, and p = 1 - (1 - /3)/n. 
Let us first restate the lower bound part of Theorem 1. 

Theorem 3.6. If/3 < 1, then 

lim lim inf J„ (?„ - yn) = 1 . 

y— »oo n—yca 

Proof. It is enough to produce a suitable lower bound on the distance of the 
distribution of St from its stationary distribution, since the chain (St) is a 
projection of the chain (Xt). 

Since 9„(s) = 0{rr^), expanding tanh[/?(5 + n"')] around in /„(*) and 
using equation (2.13) shows that, for s>Q, 

E,o[5,^i \St = s]>ps-'^- 0{n-\ (3.22) 

In 

By Remark 2, if > 

E.0 \St]>p\St\-^-^- Oin^). (3.23) 

This also clearly holds for =0 or \S t\ = n'^. (In the latter case, > 
1/n.) 
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Take the initial state 5o to be = soifi); we will specify the value of 
later. Define := |S^|p~^ whence Zq = Sq = sq. Since p"^ < 2 for large n, 
from (3.23) it follows that 

Eio[Z,+i I Zf] > Zt - , 

n 

for n large enough. Since < |5f| < 1, 

EJZ,-Z,„|Z,]. ^"'""-''"'^"^">' . ^"'""'^^^°<'^"" . (3.24) 

n n 

Applying Lemma 2.8(iii) with A = {1, 2, ...,«}, we find that 

Kso[\St\] < \so\p' + cin-"\ (3.25) 

Here and below, the constants c, depend only on /3. 

Using the variance bound Var(5,) < C2n'^ (c.f. Proposition 2.7) together 
with the inequality (3.25) shows that 

E,„ [S^] = (E,,[Sr]f + Var(5,) < sy + 2cin-'%o\p' + c.n' . (3.26) 
Taking expectations in (3.24) and using (3.26) yields 

E,„[Z, - Z,+i] < - \sIp' + 2c,n'l^\so\ + c,p-' In] + 0{n\ 

Let t* = tn - an/(l - /3). Adding the increments E^ JZ^] - Eig[Z,+i] for 
f = 0, 1, the above inequality gives that 

^0 - E JZ,] < + + ^3^^ + 0(t*n-^). 

n{l - p) n^'^ n'^{l - p) 

Since p~'* < n'^^, we deduce that 

17 r7 1 ^ ^0 ^ 2c2 log(n) o-vx 

So - E,o[Z;*] < - — - + — — +C4n ' . (3.27) 

1 — p n''-^ 

If So < (1 - /3)/3 and n is large enough, then the right-hand side of (3.27) is 
less than so/2. Thus 



By Proposition 2.7, max{Vari„(5,), Varp(5)} < Cj/n. Thus 



B/2 < E,^[Sr*] - VVar,s,(5,0, 



5/2>E^[5] + ^V^^^. 

4C5 ' 



LetTTs be the stationary distribution of (5,), and let A := [-5/2,5/2]. Then 
||P,,(5,* 6 •) - ^5IItv > ns{A) - P,„(|5,*| 6 A) > 1 - 32cle'^Vsl, 
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where the last inequality follows from application of Chebyshev's inequal- 
ity. The right-hand side clearly tends to 1 as a ^ oo. ■ 

4. Critical Case 

In this section, we analyze the mixing time of the Glauber dynamics in 
the critical case jS = 1, proving Theorem 2. We consider the upper and 
lower bounds separately. 

4.1. Upper bound. 

Theorem 4.1. ///3 = 1, then ^^ix = 0(n^'^). 

Recall the definition of tq in (2.19): tq := min{? > : |S,| < l/n}. 

Proof. We show that we can couple Glauber dynamics so that the magneti- 
zations agree in order n^^^ steps, and then appeal to Lemma 2.9 to show the 
configurations can be made to agree in another order n log n steps. 

Step 1: Our first goal is to prove that limc^oo Pcr(To > cn^^^) = 0, uniformly 
in n. 

Recall the inequality (2.17): For > n~\ 

E^[\S,^^\\S,] < |l - + ^tanh(|5,|). 

Multiply both sides above by 1{to > t} and use the fact that tanh(O) = to 
find that 

E^[|5,+i|1{to > < |i - ^jl^r|l{To > t} + ^tanh(|5,|l{ro > t}). 

Since 1{tq > f + 1} < 1{tq > t}, 

[|5,+i|l{ro > ? + 1} I 5,1 < (l - - 1 |5,|1{to >t} + - tanh(|5,|l{To > t}). 

Define ^'^ := Eo-[|5,|1{to > t}]. Take expectation above and apply Jensen's 
inequality to the concave function tanh restricted to the non-negative axis, 
to see that 

^;,i<(l--)^; + -tanh(0- (4.1) 
Thus, there exists a constant Ce > such that, if > s, then 

We conclude that there exists a time t^, = t-^(n) = 0{n) such that ^'^ < 1/4 
foralU > t^,. 
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Expand tanh(x) in a Taylor series and use (4.1) to obtain 
for t > t^. 

This shows that, for n sufficiently large, is decreasing for t > ti,. We 
will assume from now on that n is large enough for this to hold. Given a 
decreasing sequence of numbers 

\IA>bi>b2> ■■■>Q>, 

let Ui := min{? > ti, : < bi}. Since is decreasing, bi+\ < ^'^ < bi for all 
times Ui < t < Ui+i. Let bi = (1 /4)2"'. For t e (ui, 

It follows that 

16n r _^ 1 1 

M,+i -M; < —[l + 0(bfn-')\ 

Let i'o = min{/ : b, < n"'^}, where a is a parameter to be chosen below. 
If or > 2/3, then bi > n'^^^'^'^ for / < iq, for some S > 0. In particular, 
bj^ < n^'^ and 0(bj^n'^) = o{n) for / < iq. Thus for n large enough, for 
< / < io, 

32n 

Ui+i - Ui< —. 

Summing the above, 

'0-1 TO 

Z5Zn Con , 

i=Q ^io-l 

SO 



Ui, < 0{n'-^") + 0{n), 

where the second inequality follows since uq = t* = 0(n). To summarize, 
provided 1 > a > 2/3, there is a constant ci such that ^/^ < n"~^ for 
t > ciir''^". In particular, letting r„ = c\rv'~'^", there is a constant C2 > 
such that 

E^[|S;|l{To>r„}]<C2n"-i. (4.2) 
By the Markov property and Lemma 2.5, for some constant C3, 



P^(to > r„ + yn^'^ | X,-„) < 
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Multiplying both sides by 1{to > r„}, taking expectation, and then using 
(4.2) shows that 

P.(ro > rn + yn^") = 0{y~"^). 
Choosing or = 3/4 > 2/3, we see that 

P^(To > (ci + y)n^'^) = 0{y-"\ 

Step 2: Construction of coupling. We now describe how to build a Markov- 
ian coupling (X,, X,) of the Glauber dynamics such that the following holds: 
There are constants c\ > and b < \ such that, if Tmag is as defined in (3.2), 
then for any two configurations cr and a, 

P^,*(Tn,ag > Cm^''-) < b. (4.3) 

This is sufficient, since we only desire to prove f^ix = 0{n^^^). 

Fix two configurations cr and &, and suppose without loss of generality 
that \S{a)\ > |5(cr)|. Define the stopping time Tabs to be the first time the 
two chains cross over one another, i.e. 

Tabs :=minU >0 : < 

and let Gi := = |5r^|^^+i|} be the event that the two chains meet one 

step after Tabs- There is a constant C4 > 0, not depending on n, such that 

P^,^(Gi)>C4. 

On Gj, couple the two chains independently. On Gi, we divide into two 
cases: 

Case Sr„^,^\ = S T„h^+i. If this situation occurs, then couple such that the 
magnetizations continue to agree. To do so, if a site / is selected to update 
Xt with a spin S, then pick a site in X, at random from those with the same 
spin as Xt(I), and update this site also with spin S. 

Case Sr^fcj+i = -Sr„i„+h In this case, we use the reflection coupling: Sup- 
pose state / is selected to update X,, and the spin used to update is S. Then 
pick a site in % at random from those with spin -Xt{I), and update with 
spin -S. In this case, the process {St) and {St) will be reflections of one 
another for t > Tabs- 

If n is even, in either situation the magnetizations agree at time Tq, so 
Tmag ^ "^0- For ^vcn n, run the chains together after tq. If n is odd, at time 
To run the chains independently of one another for a single step. 

By Step 1 of the proof, there exists a constants c* and C6 > such that, 
for all (T, 

P<.(To + 1 < c^n^l^) > C6. (4.4) 
LetGa = {tq + 1 < c^n^'^}. 
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Let G3 be the event that the two chains couple at time tq + 1. There exists 
some C5 > not depending on n such that Va-iGs \ Gi n G2) > C5. (If n is 
even, this probability is one.) 

Then 

P<,(Gi n G2 n G3) < P^(T, < 

The probability on the left is uniformly bounded away from zero, complet- 
ing the proof. ■ 

4.2. Lower bound. 

Theorem 4.2. Suppose = 1. There is a constant Ci > such that t^^a ^ 
C,n"\ 

Proof. It will suffice to prove a lower bound on the mixing time of the mag- 
netization chain (St). 

As usual, S denotes the normalized magnetization in equilibrium. The 
sequence n^^^S converges to a non-trivial limit law as n 00. (This is 
proved in Simon and Griffiths (1973); see also Ellis (1985, Theorem V.9.5).) 
Take A > such that 

lu(\S\<An-^'^)> 3/4. (4.5) 

Take sq = 2An~^^^. Let (St) be a chain with the same transition probabil- 
ities as (St), except at so. At so, the 5 -chain remains at so with probability 
equal to the probability that the S -chain either moves up or remains in place 
at Sq. The two chains can be coupled so that St < St when both are started 
from Sq. In particular, for all s, the inequahty Fsoi^t < s) < P.s(,(S, < s) 
holds. 

Let Zt = §0 - StAT, where r := min{? > : St< An-^'"^]. Note that (Z,) 
is non-negative. 

We will now show that if Tt is the sigma-algebra generated by Zi, . . . , Zf, 
then there is a constant ca so that 

E,^[Zl,-Z^\rt]<^. (4.6) 

The equation (4.6) is clearly satisfied when Z, = 0. On the event S t = s, 
where An'^l'^ < s < sq, the conditional distribution of St+i is the same as 
the conditional distribution of given S t = s. Thus 

s^ 

E,o[5,+i I 5, = ^] = E,„[5,+i I 5, = ^] > ^ - Co-, (4.7) 

n 

for a constant cq. The inequality is obtained by expanding tanh in (2.13). 
From (4.7), it follows that 

E.„[Z,,iI^]<Z^ + -5,^ (4.8) 
n 
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We decompose the conditional second moment of Z,+i as 

E.vo[Zf^i I TA = Var(Z,+i | Tt) + (E,„[Z,+i | Tt]f . (4.9) 
Since |Z,+i - Z,\ < 2/n, 

Var(Z,^i I Tr) = Var(Z,^i - Z, + Z, | Tr) = Y^iZ^^i - Z, | Tr) < \. (4.10) 
By (4.8), for ? < T, there is a constant C\ (depending on A) so that 



E;[Z,+i I TA < Zt + 1-Z,S', + ^ < Z/ + c,n-\ (4.11) 

Using the bounds (4.10) and (4.11) in (4.9) establishes (4.6). We conclude 
that 

E.„[Zf] < CAuh. (4.12) 

Note that 

E,„[Z2] > E^JZ^Hr < ?}] > ^^.,{r < t), 
which together with (4.12) shows that 

CAt 



I (X-i / -72 , o^-O^ p;3 , t-O^f / ^2 , „ .-2 



P.o(T < < 



Taking r = (A^/4cA)n''^^ above shows that 

This, together with the bound (4.5), proves that dic^n^^^) > 1/2, where 
C3 = A^/4ca. That is, W > c^n^'^. ■ 

5. Truncated Dynamics for Low Temperature 

We now consider the case /3 > 1 . As stated in the introduction, the mixing 
time for the full Glauber dynamics is exponential in n. This is proved via 
an upper bound on the Cheeger constant, defined as 

H{x)P{x,y) 

O := mm ^ , 

A:/i(A)<l/2 //(A) 

where P is the transition matrix for the Glauber dynamics. By taking A = 
{cr : h{(t) > 0} and estimating [Z.veA,x?AM^)^(^> j)] /;"(^). whence > 1 
there are positive constants c\ and such that O < cie"'^^". The spectral 
gap of P is bounded below by C3/O (see, for example, Sinclair (1993).) The 
mixing time, in turn, is bounded below by the spectral gap (see, for example, 
Aldous and Fill (in progress).) The details of this standard argument can be 
found in the forthcoming book Levin, Peres, and Wilmer (2007). That the 
Glauber dynamics is slow mixing for yS > 1 was understood as far back as 
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Griffiths, Weng, and Langer (1966), although they lacked the tool of the 
Cheeger inequality to make a complete proof. 

Here we study the Glauber dynamics confined to the configurations where 
the magnetization is non-negative, and show that the restricted Glauber dy- 
namics has a mixing time of order n log n. 

We remind the reader of the exact mechanism for restricting the dynam- 
ics. The usual dynamics are run from a state with non-negative magnetiza- 
tion. If a move to a state r] is proposed, and rj has negative magnetization, 
then the chain moves to -77 instead. 

To establish an 0{n log n) upper bound on the mixing time, we need to 
estimate the hitting times of the normalized magnetization chain. 

Lemma 5.1. Let /3 > 1. Let s* denote the unique positive solution to 
ianhifis) = s, and for a > define 

T* = r*ia) := mf{t >0:Sl<s* + an^'^]. (5.1) 

There exists a constant c > 0, depending on a and yS, such that 

lim Po-(t* > cn log n) = 0. 



(s - s*). (5.2) 



Proof. Let 7* := yScosh ^(fis*). First, we show that 

EASt,,-s*\St = s]< 1-^i::^ 

L " 

By Remark 3 and (2.17), for S+ > 1 /n 

\st^, - 5; 1 5;1 < - [tanhOes;) -S^]. 

Since > 1, it follows that y* = /3 cosh~^ (J3s*) < 1. By the mean- value 
theorem, for y > 0, 

tmh[fi(s* + y)] - tanhO05*) = — — — y, 

cosh (s) 

for some s e [s* , s* + y]. Since cosh(x) is increasing for .x; > 0, the right- 
hand side is bounded above by y*y. Thus, for y >0, 

tanh[J3(s* + y)] < s* + y*y. (5.3) 

Hence, 

EASU - 5; I 5; = ^] < -(^ - s*) ^^~^*\ 

n 

from which (5.2) follows. 
By (5.2), 

(1-7*) 



Yr-.-- 



1 - 



-t 



(St-s*) 
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defines a non-negative supermartingale for f < t*. By optional stopping, 

1 > EAYr*^] > [(1 - (1 - r*)/n)-'^^*(5,^.^, - s*)] 

>Cin-"\l-{l-y*)/nrF^(T* >t). 
Hence Po-(t* > t) < Cin~^^\l - (1 - y)/n]', and the lemma is proved. ■ 
Proposition 5.2. Let /3 > 1. For > 0, if 

= T^(C3) := min{? >0 : St>s* + c^n^'\ 

then 

Eo[tJ = 0(nlogn). (5.4) 

Proposition 5.2 is proved in Section 5.2. Meanwhile, we state and prove 
Theorem 5.3 below, which establishes the upper bound. 

Theorem 5.3. Let^ > 1. There is a constant cifi) so that t^a{n) < c(J3)n\ogn 
for the Glauber dynamics restricted to Q.^. 

Proof. We show that there is a coupling (X/^,^/^) of the restricted Glauber 
dynamics started from states a and & such that, if r^ag is the first time t 
with = 5/", then 

lim sup Po-,a-(Tmag > cn log n) ^ as c ^ oo. 

n— »oo 

An application of Lemma 2.9 will then complete the proof. 

By monotonicity, it is enough to consider the the starting positions and 
1. The "top" chain with starting position 1 we denote by (5f), and the 
"bottom" chain with starting position we denote by (5f). Let /i^ be the 
stationary distribution of the restricted magnetization chain, and let (5,) be 
a stationary copy of the restricted magnetization chain, that is, started with 
initial distribution ji^ . 

Initially, all the chains are independent of one another. Given constants 
ci < C2, let 

Ti = \mn{t > : < s* + Cin'^''-], 

T2 = min{f > : 5f > + C2n^'^]. 

Suppose that ti < T2. On the event 5,-, > s* + C\n~^l^, for ? > ti we couple 
together monotonically the S -chain and the S ^-chain (that is, such that S , > 
S] for all t >Ti), and continue to evolve the 5^-chain independently of St 
and 5f . On the event < s* + cin"^^^, we continue to run all three chains 
independently. Then at time T2, on the event that < i'* + C2n~^^^, couple 
together all three chains monotonically (so that 5 f < S t < 5f for all t > T2). 
If 5^2 > ■^^ + '^2, just let the chains run independently. The case T2 < Ti is 
handled analogously. 
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Note that, since (S t) is independent of (S J) until after time ti , the random 
variable Sn is independent of ti and hence still stationary. 
Let C3 > be a constant, and define events Hi , H2 by 

Hi = {ti < csnlogn} n > ^* + Citi'^^^}, 
H2 = {t2 < c^nlogn} n {Sr, < s* + C2n-^^^.}. 

Then 

P^AHd < P^Ari > c^n log n) + /z+(0, s* + Cin-'l\ (5.5) 

and 

V^^Hl) < P,,^(T2 > c^nXogn) + n\s* + C2n-'l\ 1). (5.6) 

Now observe that on the event Hi C\H2 the chains (5 f ) and (5 f ) have crossed 
over by the time CT,n log n, and that by (5.5) and (5.6), 

P^,a-(//i n //2) > 1 - Pcr,^(ri > csnlogn) - P^,5-(t2 > c^nlogn)- ix^{r), 

where / = (5* + CinT^^^, s* + C2n~'^^). 

Since, as a consequence of Theorem 2.4 of Ellis, Newman, and Rosen (1980), 
the stationary magnetization satisfies a central limit theorem, ^x^{F) < 1 
uniformly in n. Further, 

lim Po-,a-(Ti > c^n log n) = and lim Po-,a-(T2 > c^n log n) = 0, 

fl— (OO ' (I— »oo 

by Lemma 5.1 and Proposition 5.2, respectively. Hence the probability that 
S ^ and S ^ will have crossed by the time c^n log n stays bounded away from 
as n ^ 00. 

Finally, observe that, whenever the two chains cross, they coalesce with 
probability bounded away from uniformly in n, which completes the 
proof. ■ 

5.1. Hitting times for birth-and-death chains. A birth-and-death chain 
on {0, 1, ... ,A^} is a Markov chain (Z,) on with transitions Z,+i - Z, con- 
tained in the set {-1,0, 1}. 

This section contains a few standard results concerning the hitting times 
of birth-and-death chains. We shall use these in the proof of Proposition 5.2 
in the next section. 

Define 



Pk 


= P(Zr+i 


-z, 


= +\\Z, = k) 


k = 0,\,- 


..,N 


-1, 


qk 


= P(Zr+i 




= -l\Z, = k) 


k = 


1,... 


■,N, 


rk 


= P(Z,+i 


-z, 


= 0\Zt = k) 


k = 


0,... 


,N. 
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Clearly, pt + qk + ft = ^ for all k if we define qo = Pn = 0. Using n to 
denote the stationary distribution of the chain, we have 



^(1) = c, 

k 



k 

7r(fc) = Cp,,,,n^, k=l,...,N, 



where Cp^^j- = [1 + 2Li Pj-il/^ ^ is a normalizing constant. 

Now, let £ < N he a positive integer, and let Z|^^ be a restriction of Z, 
to the set {0, . . . ,£}. In other words, when at k e {0,. . . ,£ - 1}, the chain 
makes transitions from k as the original chain, but when at it moves to 
i - I with probability qt and stays at £ with probability pi + rg. Let n^'^^ be 
the stationary measure of Z^f\ It is easy to verify that there is a constant 
Cp such that 

n^^\k) = Cl^^Mk) fork = 0,1,... J. 

In other words, under the stationary measure of the restricted chain, the 
states 0, 1, . . . , ^ each have the same relative weights as in the unrestricted 
chain. 

For A: G {0,1,..., A^} let 

Tk = mf{t >0:Z, = k}, 
tI = mf{t >0 : Zt = k}. 

Then (see for instance Levin, Peres, and Wilmer (2007)) for = 0, 1 , . . . , N- 
1, 

= Ef[T^] = 1 + ^,E,_i(T,). (5.7) 

In the above, Ey and respectively denote the expectation operators corre- 
sponding to the unrestricted and restricted chain starting in j. We shall now 
apply identity (5.7) to the Glauber dynamics magnetization chain. 



5.2. Hitting time for magnetization. 

Proof of Proposition 5.2. Here it is more convenient to work with M, = 
nS{X'l)l2, which is a birth- and-death chain with values in {0, ...,n/2 - 
l,n/2}. Note that, if n is odd, this chain is not integer-valued, but this 
causes no difficulties, as one can simply shift all states by -1/2. 

Let £* = lns*\. Let c > be a constant. Also, throughout the cal- 
culation, C will denote a generic positive constant whose value may be 
adjusted between inequalities. In the notation of Section 5.1, we have for 
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The probability of moving left, qe, is bounded away from 0, uniformly in 
£ 6 {1, . . . ,n/2}. Consequently, writing € = nx and j = ny, we obtain the 
upper bound 

(„,:„,,) exp(/i2„/) 
i^t-ilTe] < C — . 

Applying Stirling's formula, the right-hand side is bounded above by 
Z^=o(l +y)-^'^^'^"'\l-2y)-('-^y^"/\l-4y^)-"^exp{2j3ny^) 
^ (1 + 2x)-(i+2;c)«/2(l - 2x)-(i-2^)"/2(l - 4x2)-i/2 exp (2finx^) ' 
which can be rewritten as 

Z-=oexp[-n/Cy)](l-4/)-i/^ 
exp (1 - 4x2)-i/2 



= cJ^exp[n(/W-/0;)]|3-^j 



where 

/(z) = i(l + 2z) log(l + 2z) + ^(1 - 2z) log(l - 2z) - 2fiz\ 

Since £/n < (£* + 0( y[n))ln < 1 uniformly in n, we can bound 

l-4x^^^^^ 



sup sup - — -^1 <C. 

n 0<y<s*=f*/n \ t " / 

It follows that the behavior of each term in the sum is dominated by the 

behavior of the exponential factor exp [n(f(x) - f(y))], and so it is enough 
to upper bound the expression 

J]cxp[n(f(x)-f(y))]. 

j=o 

We then need to look for stationary points of / in the interval [0, 1] ; we have 
/(z) = log(l + 2z) - log(l - 2z) - 4;Sz 



so /'(z) = if and only if 
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or, equivalently, 

2z = tanh(2/3z). 

When P < \, the unique maximum of / is at jc = 0. When jS > 1, there is 
a local maximum of / at 5 = 0, and as mentioned earlier, there is a unique 
< 5* < 1 minimizing /. As before, we write i* = \_ns*}. 
By the above, when x < s*, 

t 

E,-i [T,] < C J] exp - f{y))\ , 

j=o 

and f(x) < fiy) for all y < x. 

Throughout the calculation below, we shall use the fact that /'(j) < for 
all y € [0, s*), and that the second derivative f"(y) exists and is uniformly 
bounded in that range, as i'* < 1/2. 

Suppose X = 0{rr^l^), i.e. € = 0( y/n). Then 

£ 

E^-i[T,] < C 2 exp [2f(x)(nx - ny) + 0(n(x - yf)] 

7=0 

£ 

<cj]cxp[(f'(e/n)(£-j)] 

7=0 

valid for 1 < £ < Ci yjn. The final bound is valid as f'(£/n) < 0, and so 
each term is bounded by a constant. 

Similarly (taking Ci = 20) we have, for 20 V" < ^ < ^*/2, 

£ 

E,-i [Td < C J] exp [nce,y)(€ - j) + 0{n{x - yf)] 

7=0 

£ 

<cj;exp[/'(Q,,)(^ -;•)], 

7=0 

where c^^y is between x and y (we could take q j, = x, for each by the 
uniform boundedness of the second derivative). There exists a constant 
Ci > such that, if j > ijl, then /'(q^) < -citjn. Then there exists a 
constant C2 > such that, for j < £12, 



f{jln)-f{£ln)<-C2. 
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This in turn implies that the sum of remaining terms is negligible. More 
precisely, 

^ exp [n(f{£/n) - fij/n))] < nexp(-C2n). 

;=o 

It follows that 

e 

Ef-iM < ^ exp[-ci^n~V^- j)] +nexp(-C2n) 

;=L^/2J 

- 7 77^ + nexp(-C2«) 

1 - exp(-cii/n) 

Cn 

for some constant C > 0, uniformly in n. 

Now suppose that t* jl < I < - 20 ^Jn. Then, for some constant 
c\ > 0, f'{ce^y) < -Ci{€* - £)/n, as long as j = yn > £/2. Also, there exists 
a constant C2 > such that, for j < £/2, 

fU/n)-f{£/n)<-C2, 

and so the contribution due to the terms with j < £/2is negligible. 

Then a calculation similar to that for 20^[n < £ < t* 12 above implies that 

there is a constant C > such that 

„ Cn 
E,-,[r,] < — , 

uniformly in n. Similarly, \i €* -20^[n < £ < \ns* + c ^/n], then we see that 

E,_i[t,] = O(V^). 

Summing over £, we obtain an upper bound on the expected hitting time 
of [ns* + c V«l starting from 0, as follows: 

[=0 

( n t*ll 



< C 



< C(n + nlogn), 

where C is once again a generic constant, and was changed to 2C in the last 
inequality. ■ 

Related results on the magnetization chain can be found in Olivieri and 
Vares (2005). 
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5.3. Lower bound. 

Theorem 5.4. Assume that /3 > 1. For the Glauber dynamics restricted to 
configurations with non-negative magnetization, t^i^(n) > (l/4)nlogn. 

The Glauber dynamics restricted to configurations with non-negative mag- 
netization will be denoted by 

Proof. Recall again that s* is the unique positive solution to tanh(J3s*) = s*. 

Since we are proving a lower bound, it suffices to consider any specific 
starting state; we take to be the all plus configuration. 

We let (X:I',X^) be the monotone coupling, where X^ is the all plus con- 
figuration and Xq has the stationary distribution /i^. We write Pi,^+ and Ei p+ 
for the probability measure and expectation operator on the space where 
is defined. 

LetS(o-) := {i : (t(i) = -1}, and 5(cr) := \S(cr)\. 

By the central limit theorem for the stationary magnetization, (c.f. Ellis, 
Newman, and Rosen (1980)), for some < Ci < 1, 

Pi„. (5(1^) < cm) = lu-'iicr : B(cr) < c,n}) = o(l). 

Let Nf be the number of the sites in S(X^) which have not been updated 
by time t. By writing A^, as a sum of indicators, 

Ei,^. [n, I B(X^)] = B{X^)[l - n-'y, 

and so, for some C2 > 0, 

Ei,^. [a^,* I 5(1^)] > C2B{X^)n-'/\ 

where t* = (l/4)nlogn. Also, since these indicators are negatively corre- 
lated, Varip+(//f) < n for all t. Applying Chebyshev's inequality shows that, 
for some C3 > 0, on the event {B(X^) > cin), 

Pi./.^ (Nr* <c,n'"\B(X;)) = o(l\ 
where the o(l) bound is uniform in B. We conclude that 
Pi,^. (Nr; < c,n'") < Pi,^. (B(X^) < c,n) 

+ Pi,^+ (a^,* < C3n^'^ and 5(1+) > c^n) 
= o(l). 

Suppose now that A^,* > c^n^^'^. It follows that St* > St* + c^n'^^^ for 
some C4 > 0. Thus, if St* < s* + c^n'^l^ for a small constant cj > 0, then 
St* < s* + (C5 - c^)n'^'^. Therefore, 

Pi,^. (St* <s* + c,n'l') < o(l) + Pi,^. (a^,* > c^n^l' and St* < s* + Csn'") 

< o{\) + Pi,^. (St* <s*+ (C5 - c,)n-'l') . 
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Again by the central limit theorem, the probability on the right-hand side 
above tends to as n ^ oo, provided we choose C5 < C4. 

On the other hand, appealing one final time to the central limit theorem, 

: S(o-)>s* + C5n-''^}) = o(l). 

Consequently, 

dniO>Pi,,^(s,*>s*+csn-"^) 

-lu^({cT : 5(cr)>5*+C5«-i/^}) 
= l-o(l), 

and so ?mix(«) > (l/4)n log n for n large. ■ 

6. Conjectures 

We believe the results proven in this paper should be generic for Glauber 
dynamics on transitive graphs. 

To be concrete, consider the J-dimensional torus (Z/nZY. Let /3c be the 
critical temperature for uniqueness of Gibbs measures on Z"^. 

We make the following conjectures: 

(i) For yS < Pc, there is a cut-ofF. 

(ii) For /3 = fie, the mixing time is polynomial in n. A stronger conjecture 
is that there is a critical dimension dc such that for d > dc, the mixing 
time is 0(\V„\^'^). 

(iii) For > fie-, if the dynamics are suitably truncated, the mixing time 
is polynomial in n. A stronger version is that again there is a critical 
dimension dc such that for d > dc, the mixing time is 0(\V„\ log \Vn\). 
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